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INTRODUCTION 

Physical  vapor  deposition  (PVD)  is 
currently  being  explored  as  an  alternative  to 
electrodeposition  for  coating  the  bore  of  large 
caliber  cannon.  PVD  is  an  efficient  and 
environmentally  friendly  means  of  producing 
protective  coatings.  A  number  of  experimental 
PVD  systems  have  been  developed  using 
cylindrical  magnetron  sputtering  (CMS)  to 
provide  a  thick  metal  coating  on  the  bore  of 
cannon.  However,  adherent,  erosion  and 
corrosion  resistant  coatings  are  critical  to  the 
perfonnance  of  many  weapon  systems  and  the 
process  control  parameters  that  produce  the 
optimal  coating  properties  have  not  yet  been 
identified.  A  number  of  promising  coatings  have 
been  produced,  but  a  significant  number  of 
costly  experiments  are  required  before  these 
parameters  can  be  established.  Therefore,  there 
is  a  need  for  further  guidance  in  selecting  the 
experiments.  A  model  capable  of  predicting 
important  coating  properties  such  as  adhesion, 
cohesion,  density,  compositional  variation,  and 
uniformity  is  desirable,  but  currently  not 
available. 

Figure  1  shows  the  interaction  of  the 
components  envisioned  for  a  comprehensive 
model  of  the  CMS  system.  The  DC-discharge  is 
simulated  using  object-oriented  particle-in-cell 
techniques  (OOPIC)  [1]  to  obtain  the  sputtered 
particle  kinetics  and  flux  distribution.  The  CMS 
magnetic  field  distribution  is  obtained  from  a 
finite  element  magnetic  field  model,  FEMLAB 
[2].  The  plasma  simulation  provides  the 
sputtered  particle  kinetics  that  is  used  by  a 
molecular  dynamics  simulation  (XMD)  [3]  to 
obtain  the  cohesive  and  adhesive  properties  and 
phase  of  the  coating  material.  Finally,  A  PVD 
feature  scale  model  (PVDPro)  [4]  uses  the 
particle  kinetics  and  flux  distribution  to  predict 
the  evolving  grain  structure,  unifonnity, 
compositional  variation,  and  roughness  of  the 
coating. 

In  this  investigation,  samples  were  prepared 
for  analysis  using  a  planar  target  and  substrate. 
PVDPro  models  cylindrical  sputtering  chambers 
with  planar  targets  [5]  and  substrates,  so 
FEMLAB  and  XOOPIC  simulations  were  not 


required  to  obtain  the  flux  distribution  used  in 
simulating  grain  nucleation  and  growth. 
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Figure  1.  Components  of  CMS  model 

BACKGROUND 

Physical  vapor  deposition  produces 
structures  that  are  topologically  and 
morphologically  metastable  [6].  Metal  vapors 
condense  into  a  fine-grained  crystalline  fonn  on 
a  cold  substrate  resulting  in  evolving  structures 
that  are  inherently  more  disordered  and  less 
relaxed  than  those  produced  by  other  coating 
processes.  The  size,  shape,  and  distribution  of 
these  grains  affect  the  material  properties  and 
therefore  the  performance  of  a  coating  [7].  The 
grain  structures  are  a  function  of  self¬ 
shadowing,  nucleation  phenomena,  surface 
diffusion,  ion  bombardment,  and  resputtering. 
They  are  also  a  function  of  process  control 
parameters  such  as  the  angular  distribution  of 
the  incoming  flux,  deposition  geometry,  target 
material,  and  sputter  gas  pressure. 

Microstructural  features  can  be  quantified  in 
tenns  of  self-affine  scaling  parameters  [8,9]  to 
provide  details  on  the  intrinsic  structure  of  the 
coating.  The  time  evolution  of  these  scaling 
parameters  uniquely  defines  the  dynamics  of  the 
deposition  process.  These  scaling  measures 
provide  insight  into  the  growth  processes  of 
coatings  and  a  means  of  calibrating  simulations 
using  experimental  data. 

APPROACH 

The  evolution  of  the  grain  structure  in 
the  simulations  and  experiments  was  assumed  to 
adhere  to  dynamic  scaling  theory.  This  approach 
is  used  to  characterize  the  evolution  of  rough 
surfaces  in  terms  of  a  growth  exponent  (a  or  H), 
roughness  exponent  (3,  and  dynamic  exponent 
1/z  [8,9].  Therefore,  the  coating  morphology 


was  assumed  to  be  a  homogenous,  self-affine 
structure  with  statistical  invariance  under 

anisotropic  dilations.  This  assumption  applied  to 
the  AFM  data  in  view  of  the  scan  size  used  in 
the  measurements,  and  to  the  simulated 

microstructure  in  the  model.  The  coatings  in  this 
study  are  single-valued  surfaces  in  E  ,  which  are 
statistically  invariant  under  transfonnations  of 
the  fonn: 

{x,y,z}  {Xx,Xy,?tHz}  (1) 

where  H  describes  the  anisotropic  scaling.  In  a 
more  general  case,  the  y  coordinate  would  have 
a  coefficient  A,K ,  and  in  the  case  of  self-similar 
fractals,  H  would  have  a  value  of  1.  The 
invariance  expressed  in  equation  1  implies  that 
any  point  on  a  self-affine  surface  can  be 

represented  in  the  fonn  {r,h(r)},  where  the 
height  function  h(r )  is  a  single  valued  function 
ofr=  {x,y}. 

Anisotropic  scaling  was  measured  in  terms 
of  the  roughness  exponent  H,  parallel 
correlation  length  q//  (7),  and  the  perpendicular 
correlation  length  qffi).  These  scaling 
parameters  provide  infonnation  relating  to  the 
intrinsic  properties  of  the  surfaces  over  a  range 
of  scales  after  t  minutes  of  sputtering.  The 
anisotropic  scaling  described  by  H  applies  over 
a  scaling  range  that  is  measured  in  terms  of 
c,//  (/).  The  parallel  correlation  length  is  the 
distance  beyond  which  there  is  no  correlation  in 
heights  between  points  on  the  surface.  The 
perpendicular  correlation  length,  ^±(0, 
characterizes  fluctuations  in  the  growth 
direction.  It  is  related  to  the  RMS  surface 
roughness,  a (7),  by  q±(7)=  V  2  *  a (7).  The  values 
of  H,  \u  (7),  and  q±(7)  were  determined  using  a 
generalized  form  of  the  height  correlation 
function  [9],  Ch(r,t)  =  ([ h(r0+r ,  t+t0)  -  h(r0t0)f ) 
to,ro  where  h(r,t)  is  a  single  valued  height  of  the 
surface  at  location  r  at  time  t.  H  is  detennined 
from  C/,(r,t)  assuming  that  the  surface  is 
continuous,  but  not  necessarily  differentiable 
[10].  This  implies  that 

|  h(r+8,t)  - h(r,t )  |  — >  0  as  8  — >  0,  however, 


lim  [^r  +  <^ — h(r,  fi]  may  nol  cxjs^  Therefore  a 
<?->o  S 

class  of  functions  is  introduced  to  describe  the 
surface  where: 

|  h{r+  S,t)-h(r,t)  |  SH.  The  derivative  of  h(r, t) 
behaves  as  §H_1  which  is  finite  at  H=l.  At  H  =  0, 
h(r,t)  is  no  longer  continuous.  Therefore,  the 
range  0  <  H  <  1  is  used  to  describe  a  the  degree 
of  differentiability  of  a  continuous  surface,  with 
increasing  H  corresponding  to  a  smoother,  more 
Euclidean  surface. 

Identifying  the  appropriate  scaling  region  in 
the  Ch(r,t )  data  is  critical  in  obtaining  accurate 
estimates  of  the  scaling  parameters.  The 
measured  parameters  depend  on  the  range  of 
data  used  to  fit  the  linear  region.  Therefore,  we 
developed  a  systematic  procedure  to  determines 
them  by  fitting  linear  and  polynomial  splines  to 
the  height  correlation  results  using  the  data  that 
minimizes  the  residuals  in  the  fit  [11],  We 
developed  this  technique  since  blind  regression 
fits  often  result  in  incorrect  values  [8]. 

Correlations  between  points  are  time- 
dependent  and  generally  increase  with  sputter 
time.  Increasing  grain  size  may  be  characterized 
by  growth  in  the  lateral  (^//(7))  direction  by  1/z. 
Corresponding  changes  in  the  growth  direction 
are  described  by  (3.  In  general,  due  to  anisotropy 
of  the  deposition  process,  (3  ^  1/z.  According  to 
dynamic  scaling  theory,  at  small  sputter  times, 
£//  (0,  and  q±(7)  are  given  by: 


4//  “  t 1,2 

(2) 

4  ~  / /; 

(3) 

The  scaling  exponents  H,  (3,  and  1/z  define  a 
unique  universality  class  that  is  independent  of 
the  details  of  the  deposition  process.  They 
provide  a  unique  metric  for  describing  the 
evolving  surface  structure  of  the  coating  and  for 
validating  the  integrity  of  the  simulation. 

A  planar  magnetron  sputtering  chamber 
was  employed  to  sputter  deposit  niobium  on  2 
inch  silicon  wafers  for  15,  30,  45,  and  75 
minutes.  The  material  was  sputtered  with  a  200 
Vdc  bias  voltage  at  10  mTorr  of  argon  gas 
pressure.  A  Digital  Instruments  Dimension  3100 
Series  Scanning  Probe  Microscope  [12]  was 


used  to  map  surface  structures  of  the  niobium 
coating  over  a  range  of  scales  from  10  mn  to 
5pm.  Lateral  resolution  was  enhanced  and 
lateral  forces  on  the  samples  eliminated  by 
oscillating  the  cantilever  at  its  first  bending 
mode  resonant  frequency  using  a  piezoelectric 
crystal.  In  this  "Tapping  Mode"  of  operation, 
the  cantilever  tip  lightly  taps  the  sample  during 
the  scan,  and  contacts  the  surface  at  the  bottom 
of  its  swing  [13].  The  resonant  frequency  of  the 
cantilever  was  287  KHz  and  the  probe  tip  was 
etched  silicon  having  a  radius  of  curvature  < 
lOnm  and  a  sidewall  angle  of  17°.  The  AFM 
images  in  this  investigation  are  5  pin  square 
with  a  10  mn  horizontal  resolution  and  a  0.1  mn 
vertical  resolution. 

PVDPro  employs  Monte  Carlo  methods  to 
simulate  the  sputtering  of  niobium  off  the  target 
surface  based  on  the  process  control  parameters 
selected  to  deposit  the  niobium  samples.  Models 
of  the  distribution  of  sputtered  material  were 
based  on  erosion  profile  measurements  of  a 
depleted  target  obtained  from  the  experimental 
system.  Gas  phase  collision  dynamics  were 
computed  with  Ml  forward  scattering  [5]  in 
simulating  the  transport  of  the  sputtered  material 
to  the  substrate.  The  distribution  of  energies  and 
angles  of  the  arriving  flux  was  then  used  by  a 
feature  scale  model  to  simulate  nucleation,  self¬ 
shadowing,  and  surface  diffusion  on  the  film 
and  substrate.  Particles  that  adsorb  on  the 
surface  migrate  over  a  given  diffusion  length 
before  being  incorporated  into  the  film.  The 
final  position  of  the  particle  is  the  site  within  the 
diffusion  length,  L,  which  minimizes  the  surface 
free  energy.  The  activation  energies  for  surface 
diffusion  on  the  substrate  and  film  are  not  well 
known  for  most  materials,  therefore  the 
temperature  dependence  is  not  explicitly 
modeled. 

The  feature  scale  simulation  of  grain 
evolution  is  particle  based,  so  surface  wetting 
cannot  be  expressed  in  terms  of  a  wetting  angle 
as  in  continuum  models.  It  is,  instead,  expressed 
in  terms  of  a  percentage,  with  100% 
corresponding  to  perfect  wetting. 

The  simulation  parameters  in  this  study 
spanned  the  entire  available  range  in  PVDPro. 
The  selected  values  for  surface  wetting  were 


20%,  50%,  and  100%.  The  substrate  and  film 
diffusion  lengths  were  both  set  equal  to  0.015 
pm,  0.026  pm,  and  0.036  pin.  The  topology  of 
0.10  pm,  0.20  pm,  0.40  pm,  and  0.80  pin  thick 
coatings  were  extracted  for  analysis.  This 
corresponds  to  t  =  0.14,  0.45,  0.55,  and  0.71 
minutes  of  simulated  sputtering  for  the  process 
control  parameters  selected  in  the  model.  The 
topology  of  the  grain  structure  was  then 
extracted  and  quantified  in  terms  of  self-affine 
scaling  exponents.  The  results  were  then 
compared  to  experimental  values  to  determine 
the  optimal  parameters  to  use  in  the  simulation. 

RESULTS 

Figure  2  shows  AFM  images  of  the  niobium 
coating  after  t  =15,  30,  60,  and  75  minutes  of 
sputtering  and  figure  3  shows  the  corresponding 
Ch{r,t).  The  C),(r,t)  data  in  figure  3 

indicate  that  the  grain  structure  of  the  niobium 
coating  is  essentially  Euclidian  (0.88),  with  the 
grain  size  (2*^// )  increasing  from  188  mn  to  368 
mn.  The  values  of  grain  size  were  validated  with 
micrographs  of  the  coatings.  Figure  4  gives  the 
time  dependence  of  ^//(t)  and  q±(/),  assuming 
that  the  evolving  grain  size  and  roughness  are 
consistent  with  dynamic  scaling  theory.  The 
value  of  (3  was  determined  to  be  0.27  and  1/z  = 
0.37  using  all  of  the  data  for  the  fit.  These 
results  agree  well  with  published  results  for 
sputtered  deposited  coatings  [14]. 

Figure  5  shows  the  grain  profiles  for  a 
simulated  0.80  pm  niobium  coating  with  20% 
surface  wetting  with  L  =  0.015pm  &  0.036  pm. 
The  particles  in  PVDPro  represent  clusters  of 
atoms  with  similar  dynamics.  Particle  size  is 
user-defined,  with  higher  resolution  simulations 
quickly  exhausting  computational  resources. 
However,  the  smallest  available  particle  size, 
5.1  mn,  was  selected  based  on  the  measured 
value  of  ^±(15)  =  37  mn  for  the  niobium 
samples.  This  imposed  computational 
constraints  on  the  size  of  the  simulations  for  this 
study.  Simulations  were  limited  to  a  maximum 
coating  thickness  of  0.80  pm. 
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Figure  2.  5  pm  AFM  scans  of  niobium  after 
t  =  15,  30,  60,  and  75  minutes  of  sputtering. 
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Figure  3.  Ch(r,t)  for  sputtered  niobium  after 
t  =  15,  30,  45,  60  and  75  minutes  of  sputtering. 

In  all  of  the  simulations,  increasing  L 
increased  the  density  and  reduced  the 
complexity  of  the  grain  microstructure.  This 
smoother,  more  Euclidean  structure  is  reflected 
in  the  increasing  values  of  H  as  given  in  Table 
1.  The  table  shows  H,  grain  size  (2*^/ ),  and  £,j_ 
for  a  0.80  pm  simulated  coating  with  20%,  50%, 
and  100%  surface  wetting  and  0.015  pm,  0.026 
pm,  and  0.036  pm  diffusion  lengths.  H 
correlates  directly  with  L  and  increases  from  a 


mean  value  of  0.54  ±  0.01  to  0.64  ±  0.02.  ^_l 
increases  marginally  with  L,  from  a  mean  value 
of  0.033  ±  0.001  to  0.036  ±  0.002  pm.  The 
effect  of  surface  wetting  on  the  analysis  is  not 
apparent. 


Figure  4.  (t)  and  ^±(t)  for  sputtered  niobium 


L  =  0.015  pm 


L  =  0.036  pm 

Figure  5.  PVDPro  feature  scale  simulation  of  a 
0.80  pm  coating  on  a  2  pm  substrate  with  L  = 
0.015  and  0.036  pm. 
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Table  1.  Scaling  parameters  for  simulated  0.80 

pm  coating  with  different  wetting  %  and 
diffusion  lengths. 

The  scaling  exponents  of  the  niobium 
coating  were  used  as  a  guide  to  detennine  the 
optimal  parameters  to  use  in  the  model.  (3  and 
1/z  were  used  to  predict  the  grain  structure  at 
the  small  sputtering  times  that  constrain  the 
model.  The  grain  structure  of  the  niobium 
coating  is  essentially  Euclidean,  so  only 
simulations  with  L  =  0.036  pm  were  considered. 
The  scaling  parameters  were  determined  for 
0.10  pm,  0.20  pm,  0.4  pm,  and  0.80  pin 
coatings  using  20%  and  100%  surface  wetting. 
The  effect  of  the  finite  AFM  probe  tip  radius 
was  also  incorporated  in  the  analysis  of  the 
grain  profile.  h(r)  was  limited  by  any 
interference  detected  between  the  simulated 
surface  and  a  probe  having  a  lOnm  radius  of 
curvature  and  17°  sidewall  angle.  This  resulted 
in  an  effective  grain  structure  with  H  =  0.81  ± 
0.01,  which  is  more  consistent  with  the 
observed  data. 

The  effect  of  surface  wetting  is  reflected  in 
the  intergranular  porosity  of  the  microstructure. 
However,  the  results  of  height  correlation 
analysis  shown  in  Table  2  suggest  that  surface 
wetting  does  not  affect  the  evolution  of  surface 
topology.  Therefore,  the  results  of  the  20%  and 
100%  wetting  simulations  were  averaged  to 
compute  the  dynamic  scaling  parameters.  This 
resulted  in  a  roughness  exponent  ((3)  of  0.39  and 
dynamic  exponent  (1/z)  of  0.29  using  all  points 
in  the  fit. 


Figures  6  and  7  shows  the  simulation  results 
and  the  predicted  values  for  ^//  (t)  and  <;%/) 
based  on  analysis  of  the  AFM  data.  The  figures 
indicate  that  the  evolution  of  the  observed  data 
and  the  simulated  microstructure  are  consistent 
with  the  assumptions  of  dynamic  scaling  theory. 
The  agreement  between  scaling  exponents  is 
reasonable  given  that  only  H  was  used  to  select 
parameters  used  in  the  model. 
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Table  2.  Scaling  parameters  for  evolving 
coating  with  L  =  0.036  pm,  and  20%  and  100% 
surface  wetting 


Figure  6.  \u  (t)  for  simulated  and  experimental 
grain  profiles. 


t  (minutes) 

Figure  7.  q±(t)  for  simulated  and  experimental 
grain  profiles. 

CONCLUSIONS 

1.  The  surface  morphologies  of  the  simulated 
and  experimental  coatings  analyzed  in  this  study 
are  consistent  with  those  of  self-affine  surface 
fractals.  In  most  cases,  a  single  exponent,  H, 
characterizes  the  scaling.  The  observed  grain 
structure  of  the  observed  and  sputtered 
deposited  niobium  is  nearly  Euclidean. 

2.  The  scaling  parameter  H,  and  dynamic 
scaling  parameters  (3,  and  1/z,  provide  a  useful 
metric  to  quantify  the  evolution  of  grain 
structure,  size,  and  surface  roughness.  The 
metric  can  be  used  to  compare  simulation 
results  with  experimental  data  and  provide  a 
means  of  optimizing  simulation  parameters 
using  experimental  measurements. 

3.  Identifying  the  appropriate  scaling  region  in 
Ch{r,t )  data  is  critical  when  measuring  subtle 
changes  in  the  scaling  parameters  of  real  data. 
Fitting  linear  and  polynomial  splines  to  Ch(r,t ) 
data  provides  a  systematic  approach  for 
measuring  the  scaling  parameters. 

4.  The  simulation  of  the  evolving  grain  structure 
of  sputter  deposited  niobium  is  consistent  with 
observed  changes  in  grain  structure,  and  both 


the  model  and  experiment  are  consistent  with 
dynamic  scaling  theory. 
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